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^^ , Summary. In the hierarchical structure formation model cosmic halos are supposed 

^^ ■ to form by accretion of smaller units along anisotropic direction, defined by large- 

[■~~. ' scale filamentary structures. After the epoch of mass aggregation (which depend 

, on the cosmological model), violent relaxation processes will tend to alter the halo 

T-H ■ phase-space configuration producing quasi-spherical halos with a relatively smooth 

^ ' density profiles. 

**! . Here we attempt to investigate the relation between halos shapes, their envi- 

^ I ' ronment and their dynamical state. To this end we have run a large (L — 500 h~^ 

Mpc, Np — 512'' particles) N-body simulation of a flat low-density cold dark mat- 
ter model with a matter density .!7m = 1 — Oa ~ 0.3, Hubble constant Ho ~ 70 
km s"'^ Mpc~^ and a normalization parameter of ag = 0.9. The particle mass is 
f^ , nip > 7.7 X 10^° h~^ Mq comparable to the mass of one single galaxy. The halos are 

"•^ ' defined using a friends-of- friend algorithm with a linking length given by / = QATv 

""^ ' where v is the mean density. This linking length corresponds to an overdensity 

i^< p/pmean — 200 at the present epoch (z = 0) and the total number of halos with 

Q ■ more than 130 particles (M > 3 x 10^^ h'^ Mq) is 57524. 

• • . 1 Halo Shapes and Environment 

•i-H . 

^^ , Halo shapes are defined by diagonalizing the moments of inertia from which 

Vh ' we derive the eigenvalues and eigenvectors of the best ellipsoidal halo fit. The 

principal axes a,b,c are related to the square root of the eigenvalues such 
that a > b > c. Our results, in agreement with other studies, indicate that 
although halos are triaxial they are significantly more prolate than oblate. 
This is quantified by using the so called triaxiality index [5] defined as T = 
(a^ — 6^)/(a^ ^ c^); which has limiting values of T = 1 and T = for the 
case of a pure prolate and oblate spheroid, respectively. Our results show that 
the fraction of halos with pronounced prolatness (ie., large Ts) is significantly 
higher than that of oblate-like halos. Overall we obtain from our simulated 
halos that (T) ~ 0.73. 
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An interesting question, especially for observational studies, is whether the 
3D halo shape distribution can be inferred from the projected (2D) shapes. 
This is an inversion problem for which, under the assumption of random 
halo orientation with respect to the line of sight and purely oblate or prolate 
spheroidal halos, there is a unique relation between the projected and intrinsic 
axial ratio distributions. This is described by a set of integral equations, first 
investigated by [6] and given by [15]: 

fil) = 3 / n „2M/s>/„2 S2TT72 ' /('?) = 1 ' 



q^ Jo (1 -g2)l/2(g2_ ^2)1/2 ' -'Vy; HJ^ (^ _ ^2)1/2 (^2 _ ^2)1/2 

(1) 

with P representing the intrinsic axial ratio while No{(3) and Np{P) the in- 
trinsic prolate and oblate axial ratio distributions, respectively. The continu- 
ous function /(g) is derived from the discrete axial ratios frequency distribu- 
tion using the so-called kernel estimators (for details see [14] and references 
therein). Inverting then the above equations gives us the distribution of true 
axial ratios as a function of f{q) (eg. [4]). Nevertheless, if halos are a mixture 
of the two spheroidal populations or they have triaxial configurations there is 
no unique inversion [12]. However, all may not be lost and although the exact 
shape distribution may not be recovered accurately one could possibly infer 
whether the 3D halo shapes are predominantly more prolate or oblate-like. 
Let us see this in more detail using our simulated intrinsically triaxial ellip- 
soidal halos which are however dominated by prolate-like shapes ((T) ~ 0.73). 
The important point here is that in order for the inverted distribution to be 
physically meaningful it should be positive for all /3's. Negative values indicate 
that the assumed model for the intrinsic halo shape is unacceptable. 

In Fig.l (left panel) we present the discrete and continuous - f{q) - dis- 
tributions of the projected in 2D axial ratios for halos of two mass ranges 
(indicated in the plot) . In the middle and right panels we present the inverted 
3D axial ratio distributions (continuous lines) for the prolate and oblate mod- 
els, respectively. It is evident that the inverted oblate-model distribution has 
many negative values which is an important indication that this model is 
unacceptable, while the opposite is true for the prolate-model distribution. 
Furthermore, we plot as histograms the intrinsic axial ratio distribution of 
the "average" prolate or oblate spheroidal fit to the 3D halo. These fits are 
realized by estimating the corresponding axial ratios by Pp = (b + c)/2a and 
/3o = 2c/ (a + b). It is evident that the purely oblate model fails miserably 
to even come close to the inverted distribution while the prolate model fits 
relatively well the corresponding inverted 3D prolate-model distribution. This 
agrees with the higher prolatness of the 3D halo shapes. We therefore conclude 
that applying the previously discussed inversion method to observational data 
(eg. [13], [11]) we can infer, even in the event of triaxial ellipsoidal halo shapes, 
the dominance of prolate or oblate- like 3D shapes. 

Another interesting fact, shown in recent high- resolution simulations of the 
concordance model, is the correlation between halo mass and halo flattening. 
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Fig. 1. Left Panel: The projected axial ratio distributions for two different halo 
mass thresholds with the nonparametric kernel estimator fit (solid line). Central 
Panel: The inverted intrinsic halo axial ratio distribution (continuous line) for the 
prolate model together with distribution of "average" prolate spheroidal fit to the 
3D halos. Right Panel: the corresponding distributions for the oblate case. 



with more massive halos being flatter (eg. [1], [8], [7]). This is counter- intuitive 
in the sense that the massive halos should collapse faster than lower mass halos 
of roughly the same formation age and thus they should have more time to 
dynamically evolve and virialize, a process that should reduce their initial 
ellipsoidal configuration. However, halo formation ages vary and secondary 
infall, which if important it could affect the halo outskirts, can produce such 
elongated halo geometry. However, even in such a case the inner parts of 
the most massive halos should be virialized and thus nearly spherical, which 
however does not seem to be the case (eg. [1] their figure 7). Note that these 
results are based on dissipationless simulations while baryonic dissipation has 
been shown to affect halo shapes, producing significantly rounder halos (eg. 
[9]). Also one may expect that halos in low-density regions, where tidal effects, 
accretion and merging are minimal to be less elongated, as indeed has been 
found (eg. [2]). 



2 Halo Dynamics & Environment 

In order to define the dynamical state of a halo we use the Z\-deviation sub- 
structure statistic [3], which looks for deviations of the local velocity mean 
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Fig. 2. Normalized axial ratio and velocity dispersion distributions. Hatched his- 
tograms represent halos with significant substructure while the dashed-line his- 
tograms represent the overall halo population. 



and dispersion around their overall halo values. For each halo particle we find 
its Nnn nearest neighbors from which we calculate their local velocity mean 
and dispersion which we then compare with their overall halo values. The lo- 
cal deviations are defined by: 6f = -^ [{viocai - w)^ + {^vJocai ~ cr„)^] while 
the individual halo Z\-statistic is given by the sum A = ^5i. It has been 
found that a robust measure of the substructure index is given when using 
Nnri = 25 [10], which is the value that we use. In figure 2 we show the velocity 
dispersion (right panel) and /3p axial ratio (left panel) distributions for halos 
in the mass range 3 x 10^^ < M < 8 x 10^^ h^^ Mq. Those that have a 
substructure index Z\ > 1.5 are shown as hatched histograms while the over- 
all halo distributions as dashed-line histograms. It is evident that dynamically 
young halos (ie., those with significant substructure) are more elongated and 
have a higher velocity dispersion than the overall halo population. 

Furthermore, we investigate the correlation between halo velocity disper- 
sion and halo mass. From our halo identification procedure we expect that 
most halos should be nearly virialized and thus their velocity dispersion should 
be strongly correlated with their mass via the virial theorem. Indeed, there 
is a strong correlation between the velocity dispersion and the halo mass, 
measured by summing the member DM particles masses, with a Spearman 
correlation coefficient of _R ^^ 0.9. We find however, that this strong correla- 
tion breaks down in the vicinity of large halo hosts. In figure 3 we show with 
open points the correlation coefficient between halo a^ and Mass as a function 
of distance to their nearest massive (M > 2 x 10^^ h~^ Mq) host neighbor. 
The monotonic drop of the correlation coefficient with decreasing halo-host 
separation is evident. We also find that there is a significant although weaker 
correlation between a^ and the substructure index (Z\), which increases as the 



Halo Shapes, Dynamics and Environment 




separation (h ' Mpc) 



Fig. 3. Spearman correlation coefficient for the a^ vs Mass and a^ vs A correlations 
as a function of host-halo distance. 



halo-host distance decreases. This probably implies that halos near large hosts 
are either disrupted due to the presence of strong tidal fields or more probably 
that the excess density of sub-halos near hosts induce strong inter-halo gravi- 
tational interactions and merging. These results have important consequences 
for observational studies and put important limits on the validity of using the 
virial theorem to estimate group or cluster masses in the vicinity of massive 
clusters. 
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